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SUPERSONIC POTENTIAL FLOW 
By Herbert C. Nelson and Herbert J. Cunningham 

SUMMARY 


A Rayleigh type analysis involving chosen modes of the panel. as 
degrees of freedom is used to treat the flutter of a two-dimensional flat 
panel supported at its leading and trailing edges and subjected to a 
middle-plane tensile force. The panel has a supersonic stream passing 
over its upper surface and, still air below. The aerodynamic forces due 
to the supersonic stream are obtained from the theory for linearized two- 
dimensional vinsteady flow and the forces due to the still air are obtained 
from acoustical theory. 

In order to study the effect of Increasing the number of modes in 
the analysis, two and then fovir modes are employed. The modes used are 
the first foxur nat\aral modes of the panel in a vacuum with no tensile 
force acting. The analysis includes these variables: Mach number, struc- 

t\jral damping, tensile force, density of the still air, and edge fixity 
(clamped and pinned). For certain combinations of these variables, sta- 
bility boundaries are obtained which cam be used to determine the panel 
thickness required to prevent flutter for any panel material and aJLtitude. 

In contraist to some previous panel-flutter Investigations, the 
results of the present analysis show that sufficiently thick panels are 
flutter free for the Mach numbers treated and sviggest that this is true 
throughout the supersonic speed range. ' ' . • . • . 

• INTRODUCTION 


The flutter of thin metal plates or panels, such as compose the 
covering or skin of missiles and other craft intended for high-speed 
flight, has recently become of increased concern. Such panels may be 
initially flat or curved arid may be small or fairly large in aspect 
ratio. In addition, they may be prestressed and will probably become 
warped in flight by aerodynamic heating. If one or more of the panels 
on a particular configuration are vibrating, the basic structure 
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supporting them can usually be considered rigid. The fixity at the 
edges of the panels ranges between clamped and pinned, depending on the 
construction. Some preliminary experimentation and analytical work sug- 
gests that this type of instability is of concern only at supersonic 
speeds . 

The problem of panel flutter embraces so many possible factors as 
to discourage general treatment, and previous papers on the subject (for 
example, refs. 1 to 7) have employed various simplifying assumptions in 
order to obtain specific solutions to what might perhaps be considered 
different phases of the problem. In aill the references cited, the main 
assumption made is that a panel and the flow over it are two-dimensional. 
Other assumptions common to the reference papers are that small-deflection 
plate theory and linearized flow theory may be used. 

References 1 to 4 examine the case of a panel buckled by a constant 
shortening and held at its leading and trailing edges, with a supersonic 
stream over its upper surface and no perturbation pressures on its lower 
svirface. In reference 1 steady-state air forces and in reference 2 quasi- 
stationary air forces, which include the first order of the frequency of 
oscillation, are used. Both these references consider the dynamic sta- 
bility of the buckled panel. Reference 3 and the more exhaustive ref- 
erence 4 examine the static stability of the buckled panel and propose 
that motion (flutter) is the result when static equilibrium is not pos- 
sible. Reference 5 and a section of reference 2 treat the case of a 
flat panel pinned at its leading and trailing edges. Reference 5 uses 
exact linearized unsteady aerodynamic forces and therefore, in contrast 
to reference 2, imposes no limitations on the order of the frequency. 

In references 1, 2, and 5 a genera.li zed-coordinate approach involving 
chosen modes of the panel as degrees of freedom is employed. Reference 6, 
on- the other hand, indicates how the problem of a vibrating membrane in 
a supersonic stream can be treated by means of Laplace transforms and 
s\iggests that similar treatment can be given to the plate problem. Ref- 
erence 7 treats the case of a two-dimensional panel on many equally spaced 
simple supports with compressible air flowing over the upper surface and 
dead air below the panel, and the results indicate that the possible 
panel instabilities are divergence for svibsonic flow and flutter for 
supersonic flow. Some questionable features of the results obtained in 
references 2, 5, and 7 are examined in the section entitled "Results and 
Discussion" in the present paper. 

In the present paper, a Rayleigh type flutter analysis is developed 
by means of Galerkin's process for a two-dimensional flat panel held in 
some manner at its leading and trailing edges and acted on by a middle- 
plane or axial force (which, in the case of tension, introduces a restoring 
force similar to that for the membrane). The upper surface of the panel 
is subjected to a supersonic stream and the lower surface to an tinconfined 
mass of stationary compressible air. The normsil modes of the panel with 
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no middle -plane force acting are used as degrees of freedom in the anal- 
ysis. As in reference exact linearized unsteady aerodynamic forces 
are employed. In the reference paper the integrals yielding these forces 
are evaluated analytically. In the present paper these integrals are 
evaluated numerically. 

Numerical results are presented in order to examine some effects of 
including two and then fovir modes in the analysis and to determine effects 
of Mach number, density of the supersonic stream, panel mass and stiffness, 
edge fixity (to some extent), structural damping, axial load, and density 
of the still air below the panel. In appendix A an alternative solution 
by means of Laplace transforms is developed for the plate problem just 
described. No niamerical results are obtained by this method, however. 


SIMBOIS 


a 

Amn7®mn^*-TTin 


c 

D 

E 

f 

F 

g 

Gmn 


speed of sound in undisturbed medium 

structural and aerodynamic integrals defined after 
equation ( 15 ) 

panel chord 

local flexural stiffness^ : 

12(1 - v 2 ) 

Young’s modulus of elasticity of panel material 
tension parameter, F jc^wa^^ 

functions defined in equation (24) 

external force per unit width acting in midplane 
of panel (tensile force positive) 

structural damping coefficient 

matrix elements defined in equation (l3) 

Hankel function of second kind, of zero order, 
(notation of ref. l8) 

aerodynamic integrals defined aifter equation (23) 
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Jp(u),Yp(u) 

k 


Kn 

L(u) ,L(u) 
Lr(u) 

®A 


M 

Ni,N2 

p(x,t) 

Pn(x) 


Pu^Pl 


Poo^PQ 


Pn(x)>Pn(x) ,Pn(x) 


1 


Sr 


t 

U 


x,y,z 


Bessel functions of order p of first and second 
kind, respectively, (notation of ref. I8) 

reduced frequency, C(o/2U 

stiffness parameter (reduced first natural fre- 
quency) , caj^^2U 

eigenvalxies defined after equations (16) and (17) 
and given for first four panel modes in table I 

aerodynamic functions defined after eqviation (21) 

aerodynamic functions defined after equation (B16) 

panel mass per unit svirface area, ar 


Mach number, U/a 

coefficients in mode -shape equations (16) and (l7)j 
given for first four panel modes in table I 

net perturbation pressure, positive downward 

pressure coefficient associated with mode shape 
Zn, defined in equation (20) 

upper- and lower-surface contributions to 
pertvirbation pressure, respectively 

pressures in undisturbed supersoniq stream and 
still-air region, respectively 

components of Pn(x) defined after equation (22) 

1 p 

dynamic pressure, — pU 

coefficients of equations (28) and (BI5) tabulated 
in appendix B for first four modes of panel with 
clamped edges 

time 

velocity of supersonic stream 

coordinates defined in figure 1 
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z(x,t) 

Z(x) 

Zn(x) 

V 

P^PO 

a 

T 

(X> 

(L 

[J 

( I 

Subscripts : 

u 

I 


5 


vertical displacement of panel 
flutter mode shape 

nth natural mode shape for panel vibrating in vacuo 
panel-air mass ratio, ny^/pc 
Poisson's ratio 

densities in undistvirbed supersonic stream and 
still-air region, respectively 

density of panel material 

local thickness of panel 

distxnrbance-velocity potential 

frequency of oscillation 

frequency associated with mode shape 

frequency parameter, 

/ccn \2 

frequency ratio squared (— ) f except in flutter 

calculations where ft = (l + ig) 

square matrix 
column matrix 

upper-surface contribution 
lower-s\n:face contribution 


Primes denote differentiation with respect to the argument. 
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ANALYSIS 

Statement of Problem 


A thin, isotropic, two-dimensional plate (beam) of constant thickness, 
as shown in figure 1, is considered herein. The plate is undergoing simple 
harmonic motion and is acted on by a middle -plane or axial force F 
(tension or compression); its upper surface is subjected to a supersonic 
stream of velocity U, pressure p^o, and density p, and its lower sur- 
face is subjected to still air with pressure pq and density pq* The 

differential equation of motion for the plate may be written as 


D ^ + mA - F + (p„ - Pq) + p(x,t) = 0 




_ p B^z 




( 1 ) 


where, in the present case, the vertical displacement of the plate z(x,t) 

may be expressed as Z(x)e^^^, o) is the circvilar frequency of oscilla- 
tion, p(x,t) is the net perturbation pressure (positive downward) arising 
from the motion of the plate, mA is the plate mass per unit surface area, 

and the local flexural stiffness D is given by Et3/i 2(1 - v^) . For 
the case where p(x,t) = 0, equation (l) may be obtained, with appropriate 
changes in notation, from reference 8. 

In the remaining development the constant-pressure term p^^^, - Pq 

of equation (l) is considered to be zero. This in no way affects the 
generality of the results for the oscillating plate, since inclusion of 
the constant-pressure term as nonzero would result only in adding a 
particular solution which represents a static vertical deflection. In 
addition, the coordinate x of eqmtion (l) is divided by the plate 
chord c and henceforth is employed in this nondimens ional sense. Thus, 

equation (1) multiplied by becomes 


4 Z"" - - — Z" + p(x,t)e-^^ = 0 (2) 


where the primes denote differentiation with respect to the argument x. 

In order to obtain a specific solution of equation (2), four boundary 
conditions are required. In the present paper, the plate is considered 
to be held at its leading and trailing edges as shown in figure 1, and 
this assvimption leads to the conditions for pinned edges: 



NACA TN 5^65 


7 


Z(0) = Z(l) = Z"(0) = Z"(l) = 0 (3) 

and for clamped edges: 

Z(0) = Z(l) = Z'(0) = Z'(l) = 0 (If) 

In a later section of the analysis the boundary- value problems, as exempli 
fled by equations (2) and (5) or (2) and (4), are solved by Galerkin's 
method. Also considered in appendix A are the solutions to these prob- 
lems by means of Laplace transforms. 


Net Perturbation Pressvtre p(x,t) 

The net pressure p(x,t), as mentioned previously, arises from the 
oscillatory motion of the plate. It is this pressure which damps, or 
in the case of flutter sustains, the oscillation. In the present paper 
the pressvire on the upper surface is obtained from the theory for line- 
arized vinsteady supersonic flow and the pressure on the lower surface 
from acovistical theory. The perturbation pressure in terms of the presr 
sures p^ on the upper svirface and pj on the lower svirface is 


P(x,t) = Pu - Pz 


where 


(5) 


and 



U ^u\ 
c Sx / 


(6) 






(7) 


From reference 9 the velocity potential for the upper surface can be 
obtained in the form 





e-i^x-|)j^ 






( 8 ) , 
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where 


_ 2m^ 

( 0 = ^ 



k = 


ooc 

2U 


Based on reference 10, the velocity potential for the lower s\arface can 
be obtained, as shown in appendix B, in the form 


r 2(OHo^^^(2kM|x-||)d| 
2 J Q 


( 9 ) 


where Hb^^Ux) is the Hankel function of the second kind, of zero 
order . 


Solution by Galerkin ' s Method 

Outline of method .- The boundary-value problems considered earlier 
(eqs. (2) and (5) for the pinned-edge plate and equations (2) and (U) 
for the clamped-edge plate) are now solved by means of Galerkin 's method. 
(a detailed account of Galerkin *s method may be fovuid in ref. 11.) As 
a first step, the flutter mode shape Z(x) is approximated by a linear 
combination of the form 


N 

z(x) = c 21 anZn(x) 

n=l 


( 10 ) 


where the coeff icients a^. "may represent complex amplitudes and. where 
the functions :Zn(x) are the mode shapes for the plate • vibrating in a 
vacuum without an axial force F acting. The function is the 

fiindamental mode shape associated with the lowest natural freqiiency 
and the remaining functions ^ 2 ^ • Zjj are consecutively the 

higher modes shapes. The shapes satisfy the pinned-edge or clamped- 

edge boundary conditions (Zn replacing Z in eqs. (3) and {k)) and the 
differential eqmtion 
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4 Zn" - an^DAZn = 0 


( 11 ) 


where o^i is the frequency of oscillation for which is the mode 

shape • 


The remainder of the Galerkin process for solving the aforementioned 
hoimdary-value problems consists of determining the coefficients ajj of 
equation (lO) in the following manner: Substitute equation (lO) into 


equation {2), replace the term Zn"" by 1^ accordeince with 

cA 

eq\iation (ll), multiply by one of the mode shapes integrate the 


result from x = 0 to x = 1, and equate to zero. When n is made 
1, 2, . . . N in succession, W linear equations are obtained which 
determine the unknowns an* These equations can be written in the form 


•^11 ^12 • 

• • Gm 


^1 


0 

G21 G22 • 

• • Gpu 


^2 


0 

• • « 




^ ^ 

> 

GnI Gn2 • 

• • Gun 


ajj 


0 

_ 

__ 


V > 




The matrix elements are given by 


where 


Gmn - H lAan 


n 




- Cum 



f 




Amn = 

f ZjnZn dx 

^ 0 



®ran “ 

/ Zm'Zn' 
^ 0 


pi . 

^mn ~ 

/ ZmPn(x)^ 

d n 


0 


( 12 ) 


( 13 ) 
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and where Pn(x) is the pressure p(x,t), obtained from equations (5) 

to ( 9 ) with Z replaced by multiplied by e-i“*/pca:?. Inequa- 

tion (15) o:q and cqa are the first and the mth natural frequency, 
respectively, of the plate with no eixial force acting. 


Flutter determinant.- The flutter condition or condition of harmonic 
vibration, which is given by the nontrivial solution for the coeffi- 
cients an, is obtained from equation (12) by setting the determinant 
of the matrix G^n equal to zero. Thus the flutter condition may be 

expressed in the form 


^11 ®12 • * • ® 1W 


G21 G22 • 



( 14 ) 


Gjji Gjj 2 • • • Gun 


Remarks on alternative procedure .- The procedure from equation (U) 
to equation (15) is, in general, not the most accurate that could be 
followed for values of F other than zero. A generally more accurate 
procedure wo^lld be to use, instead of equation (U), the differential 
equation for the panel with tension: 


^ Zn - (%^niA^n " ” ^ - 0 
c^ c2 


(15) 


When eqxiation (15) is solved, subject to the appropriate boundary con- 
ditions, the frequencies cqi are found to be functions of F for both 
pinned and clamped edges, but the mode shapes Zj^ do not vary with F 
for pinned edges. The use of equation (15) rather than equation (U) 
would mean that in eqviation (13) the term fBmn woxald not appear and 
ihe frequencies and mode shapes would be those that satisfy equation (15) • 

In the present paper, equation (ll) rather than equation (15) has 
been used for the following reasons: For pinned-edge panels there is 

no difference in the mode shapes or in the final numerical flutter restilts 
for clamped-edge panels the determination of the values of Zjj and 

that satisfy equation (15) is laborious and must be carried out for every 
desired value of F. Elimination of the term f Bmn from the matrix 
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elements, through use of eqviation (15); does not compensate for the labor 
of determining the natural frequencies and mode shapes as functions of F 
The differences in final nimerical flutter results for the clamped-edge 
panel approach zero as the nimiber of modes in the analysis is increased 
and are expected to be small even when only a few modes are used. 


Eval\iation of Terms in Flutter Determinant 

Structural integrals A ^in and and frequencies Considers 

tion wi]J. now be given to the evaluation of the mode-shape integrals and 
frequencies in the elements of equation ( 15 )» The mode shapes and 

associated natural frequencies ofa obtained from equation ( 12 ) are: 

For the pinned-edge plate, 


Zn = Nq sin KnX 



(16) 


where Kn is obtained from the frequency equation 


sin = 0 


For the clamped-edge plate, 

Zn = Nq|cos K^x - cosh K^x + N2(sin K^: 


X - sinh K^x) 


> (17) 


Oil = 


where is obtained from the frequency equation 


cos Kn cosh Kn = 1 

In equations (16) and (17) the factor Nq is used to produce vinlt deflec 

tion at the center of the plate (x = O.5) for modes that are symmetric 
about the center and at the point of maximum deflection between the 
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leading edge (x = 0) and the center of the plate for modes that are 
antisymmetric about the center. The factor N2 in equation (17) is 
established by the boundary condition requiring zero deflection at x = 1 
and is expressed by 


cos Kn - cosh Kn 
N2 


sin Kn - sinh 


(18) 


The quantities ®mn required in equation (l 4 ) 

can be determined directly from equations (16) or ( 17 )- First, however, 
values must be established for Nq and Kn in the case of the pinned- 
edge plate and for Nq, N2, and Kn in the case of the clamped-edge 

plate. Table I includes values of all these qviantities for the first four 
modes of vibration. 

The values for Ajnn shown in table I are zero when m ^ n becaixse 
of the orthogonality of the mode shapes of equations (16) and (l7)» 

For the pinned-edge case the slopes Zn' of the mode shapes are also 
orthogonal and, consequently, Enm is zero when m ^ n. For the clamped- 
edge case, even though the slopes Zn' are not orthogonal, the integrand 
of is antisymmetric about x = 0.5 when m and n are not both 

even or both odd and, consequently, Bmn is zero when m ^ n except 
for Bq^, Bjq, B 24 ^ and first four modes). 

Aerodynamic integrals Cmn-- The remaining term in the elements of 
eqiiation (l 4 ) that requires evaluation is the integral 



As mentioned previously, Pn(x) is the pressvu'e p(x,t) , obtained from 
equations (5) to ( 9 ) with Z replaced by Z^, multiplied by 
The quantity Pn(x) is therefore given by 
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The quantities Jp and Yp axe the Bessel fvinctions of order p, of 
the first and second kind, respectively. 

For convenience p^(x) is considered in three parts, namely, 

Pn(^) ~ ^n(^) ^n(^) ^n(^) (22) 

where 


P„(x) = i 


f Z„(t)L(x-{)d| + ~ Zn(x) + Zn’(x) 

Jo 


Pn(x) = — r Z„( 6 )£(x-t)« 

P Jq 


Pn(x) = Zn(|)loge(kM|x - ||)d| 


Hence, equation (l 9 ) the fonn 


Cmn = ^mn + Imn + Imn 


where 


^mn 


= f 2 aPn(x)dx 


(25) 


Imn = / 

^0 


Imn = f ZjjjPn(x)dx 
^0 
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= ?(/o’ 


-ZnU) + 


^-lS(x-i)j^ 




■n + 


- - fr 

2k dx [Jq 


-z„(t) + J. 
^ 2k 


ri^x-|)j^ 


^x-l) 

M 




/u 


£a 

p 


Z^(|)Ho(2^2kM|x-||)d| 

.2^0 Jj 


( 20 ) 


where the contributions from the upper and lower surfaces of the plate 
are designated by subscripts u and I, respectively. Upon elimination 
of the derivatives in the integrands of the upper-surface contribution 
thro\;igh integration by parts, performance of the indicated differentia- 
tion of the second integral, and extraction of the singularity at | = x 
in the Hankel fxmction of the lower-surface contribution, equation (20) 
may be written in the form 


Pn(x) = ^ 
P 


Z„U)L(x-|)d{ + i ^ z„'(x) 


u 


Po 

P 


I Zn(4)L(x-|)d| + i 
0 ^0 


Zn(l)loge(^l* - 


(21) 




where 


L(u) = 






.-icm 


L(u) = i Jo(2kM|u|) + i Yn(2kM|ul) - i loge(kM|u|) 
2 2 ^ n 


L(0) = 


Z + i 

7t 2 


7 = 0.577216 (Eiiler's constant) 
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The first integral Imn represents the effect of the supersonic stream 
passing over the upper surface of the plate; the other two integrals 
represent the effect of the still air below the plate. 

Before further development of the method of the present paper for 
determining Cnm^ the aerodynamic treatments of references 1, 2, and 5, 
which deal with the pinned-edge plate, will be examined. These refer- 
ences consider only the effects of the supersonic stream, the air below 
the plate being treated, in essence, as massless; that is, Pq/p is 

taken to be zero and the integrals I mn a^id Imn are omitted. In ref- 
erences 1 and 2 the aerodynamic effects are accounted for as if the 
integral Imn has been expanded as a power series in the frequency of 
oscillation; reference 1 retains only the steady-state or zero-order 
frequency term and reference 2 adds the first-order frequency term. In 
reference 5, on the other hand, the integral 1^^ is evaluated exactly 

with regard to the frequency. This is possible because the modal func- 
tions Zn for the pinned-edge plate are sine waves (see eq. (16)) so 
that linn obtained in terms of the functions (sometimes called 

Schwarz fxinctions) 








(A = 0,1) 


(24) 


where 



CO 


b has the four values 



A similar result could be obtained for the clamped-edge plate by approxi 
mating the modal functions Zn (see eq. ( 17 )) by a finite sine series 

R 

Zn ** ^ dj. sin ntx 
r=l 
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For either pinned or clamped edges, the argvunents a and b of the 
Schwarz functions would range from large positive to large negative 

values, particularly for the higher modes, and would thus require exten- 
sive tabulation of fQ and f^. The present paper employs the exact 
expression for the pressure term Pn(x) but, because the necessary tables 
of fo and f^ are not available, for convenience, a numerical method of 
integration is used to evaluate Pn(*) and the aerodynamic integrals Imn^ 
linn, and Imn of equation (23). 


The n\jmerical method is based on the following integration rules for 
parabolic arcs: 


f ^ y(x)dx = — 5 y(xi) + 8y(x2) - yCxj) 
J XI 12 L J 


(25a) 


/ y(x)dx = 
^X2 


12 


-y(xi) + 8y(x2) + 5y(x5) 


(25b) 


where X2 = X]^ + Ax and x^ = X2 + Ax. The range of Integration in 
equation (l9), 0 = x^ 1, is, for convenience, divided into an even number 

of eqvial segments. From the standpoint of accuracy the number of segments 
needed depends on the number of nodes in the highest mode and on the value 
of ^ for which Pn(x) is evaluated. For the numerical applications of 
the present paper 10 segments were found to be adequate, and the method 
of integration is illustrated for this number of segments in the equations 
to follow. 


The use of 10 segments would, in general, require the determination 
of Pn(x) in equation (19) at 11 points on the plate. However, the 

integrand of C mn at x = 0 and x = 1 is zero since the mode shapes 
are zero at these points, and therefore Pn(x) need be evaluated at only 
the 9 interior points (equally_^spaced between x = 0 and x = l) . 

The values of the terms Pn, Pn, and Pn of equation (22) for these 

points may be arranged in matrix form eis follows: 
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Pn(*l) 


8l( 0) -L(-.l) 0 

0 0 0 

0 

0 

0 


Zn(-l) 




Zn’(.l) 

Pn(-2) 


13L(.l) 7L(0) -L(-.1) 

0 0 0 

O 

0 

0 


Zii(.2) 


Zn(-2) 


Zn'(.2) 

Pn(-5) 


15L(.2) 12L(.l) 7L(0) 

-L(-.l) 0 0 

0 

0 

0 


Zn(-3) 


Zn(-5) 


Zn'(-5) 

Pn(.M 


13L(.3) 12L(.2) 12L(.l) 

7L(0) -L(-.l) 0 

0 

0 

0 


Zn(-M 





Pn(-5) 

1^ 

13L(.1+) 12L(.3) 12L(.2) 

12L(.l) 7L(0) -L(-.l) 

0 

0 

0 


Zn(-5) 


2n(-5) 


Zn'(*5) 

Pn(-6) 

" 120p 

13L(.5) 12L(.1+) 12L(.3) 

UL(.2) 15L(.l) Ul(0) 

0 

0 

0 


Zn(-6) 


Zn{ *6) 

l+k20 

Zn'(.6) 

Pn(-7) 


13L(.6) I2L(.5) 12L(.1+) 

HL(.3) 14U-2) 12L(.1) 

5L(0) 

0 

0 


Zn(-7) 


Zn(-7) 


Zn'(.7) 

Pn(*6) 


13L(.7) 12L(.6) 12L(.5) 

Ul(.U) ll+L(.3) Hl(.2) 13L(.l) 

5L(0) 

0 


Zn(-8) 


M-8) 


Zn'(.8) 

Pn(.9) 


13L(.8) 12L(.7) 12L(.6) 

DX(.5) ll+L(.U) Ul(.3) 12L(.2) 13L(.l) 5L(0) 


Zn(-9) 


Zn(-9) 


Zn'(-9) 




L( 0 ) L(.l) L(. 2 ) L(. 3 ) L(.l^) L(- 5 ) £(- 6 ) £{.7) £(- 8 ) 


13Zn(.l) 

Pn(-2) 


L(.l) L( 0 ) L(.l) L(. 2 ) £(. 3 ) £(.M £(.5) £(. 6 ) L(.T) 


12Zn(.2) 

Pn(-3) 


L(. 2 ) L(.l) L( 0 ) L(:i) £(. 2 ) £(.3) £(. 4 ) £(.5) £(. 6 ) 


12Z„(.3) 

Pn(-‘») 

_ ^ 
” 120 

L(. 3 ) £(.2) £{.1) £(0) £(.1) £(.2) £(.3) £(- 4 ) £{.5) 


n 2 ^{. 4 ) 

Pn(*5) 

£(. 4 ) £(.3) £(.2) £(.1) £(o) £(.1).£(.2) £(.3) £(.U) 


l4Zn(.5) 

Pn(-6) 

£(.5) £(. 4 ) £(.3) £(-2) £(.i) £(o) £(.i) £(.2) £(.3) 


UZn(.6) 

Pn(.7) 


£(.6) £(.5) £(.M £(.5) £(.2) £(.i) £(o) £(.i) £(.2) 


12Zn(.7) 

Pn(«8) 


£(.7) £(. 6 ) £(.5) £(. 4 ) £(.3) £(.2) £(.i) £(o) £(.i) 

i 

12Zn(.8) 

Pn(-9) 


£(. 8 ) £(.7) £(. 6 ) £(.5) £(A) £(.3) £(.2) L(.i) L(o) 


13Zn(-9) 


Pn(-l) 


Li(.l) L 2 (. 1 ) L 3 (. 1 ) L4(-1) L 5 (. 1 ) L 6 (. 1 ) l 7 (.l) La(.l) L 9 (.l) Lio(.l) 


«l 

Pn(.2) 


Li(. 2 ) L 2 (. 2 ) L 3 (. 2 ) £j*(. 2 ) L 5 (. 2 ) L 6 (. 2 ) lrj{.2) Iq(. 2) lg(.2) Lio(-2) 



Pn(.5) 


Li(.3) L2(.3) L3(.3) L4(.3) L5(.3) ^(-3) L7(.3) L8(.3) ^(-J) Lio(.3) 



Pn(-^) 

- pq/p 

1 ” 4 

Li(.4) L2(-4) L3(.4) lj*(. 4) L5(.4) L6(-4) L7(.4) La(.4) L9(.4) Lio(.4) 


84 

Pn(-5) 

Li(.5) L 2 (. 5 ) L 3 (. 5 ) L4(.5) L 5 (. 5 ) Ui-^) L8(.5) Lg(<5) Lio(.5) 


85 

Pn(-6) 


Li(. 6 ) L 2 (. 6 ) L 3 (. 6 ) L4(.6) L 5 (. 6 ) l 6(-6) La(.6) L 9 (. 6 ) Liq(.6) 


®6 

Pn(-7> 


Li(.7) L2(.7) L3(.7) L4(.7) L5(.7) L6(-7) L7(.7) La(.7) L9(-7) Liq(.7) 


57 

Pn(-8) 


Li(. 8) L2(.B) L3(.8) L4(.8) L5(.8) L6(.8) £^(.8) La(.8) L9(.8) Lio(.8) 


«8 

Pn(-9) 


Li(.9) L2(.9) L3(.9) lk{-9) L5(.9) L6(-9) Ly(.9) La(-9) l9(.9) Lio(.9) 


69 
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where L(u) and L(u) sire defined after equation (21) and Lr(u) 
and Sj. are defined in appendix B (eqs. (B16) and (B15))- The rows of 
the sqiiare matrix in equation (26), down to the row pertaining to Pn(*5) 
were obtained by applying the integration rule given by equation (25a). 

In the remaining rows the contributions from the region x > O.5 were 
obtained by applying eqviation (25b). In equation (27) the integrating 
factors, multiplying Zn in the column matrix were obtained by using 
equation (25a) between x = 0 and x = 0-5 and equation (25b) .between 

X = 0.5 and X = 1.0. In appendix B the singular integral Pn(x) pre- 
sented after equation (22) is evaluated and the form of P^ leading to 
equation (28) is derived. 

Equations (26), (27)., and (28) are siamned in accordance with equa- 
tion (22) to obtain the column matrix |Pn(x)|- • By '*^se of this column 

matrix and the integration rtiles of equation (25) > the aerodynamic term 
C mn is obtained in the form 



l8 


NACA TN 5465 


•'mn 


1 

120 


15Pn(*l)Zm(-l) + 12pn(.2)Zm(.2) + 12pn( .5)Zm( .3) + 


l^Pn( •5)Zin( *5) + HPn( 

12Pn(-T)Zm(-T) + 12Pn(.8)Z^j^(.8) + 13Pn( .9)Zm( *9) 


(29) 


where the integrating factors 13 , 12, . . . 12, I 3 were obtained in the 
same manner as those in equation ( 27 )- By means of equation (29), C mn 
can be evaluated for a given edge fixity and for particvilar values of M, 
k, and Pq/p- 


Solution of Flutter Determinant 


As previously stated, the conditions for flutter are determined 
from the nontrivial solutions of equation (l4). Since equation (l4) is 
complex, it may be solved directly for one complex imknown or two real 
unknowns. For a specific edge fixity the variables (see eq. (15)) in 
equation (l4) are l/p (the inverse of p is preferred because p 
becomes infinite for p = 0) , D, f, M, k, and Pq^/p- In the present 

paper it is convenient to interpret the D of equation ( 15 ) as the complex 
quantity f— ) (l + ig) rather than f— ^ , where g may be regarded 

\(o/ \o)y 

as a structural damping coefficient. (For this use of g, see, for 
example, refs. 12 and 15 .) Each of the various quantities on which equa- 
tion (l4) is dependent was varied to some extent, as will be discussed 
in the next section. A particular calculation was performed by setting 
values for l/p, M, k, f, and Pq/p nnd solving for D. Then, because 

it was one of the more easily varied parameters, l/p was changed and 
again fl was solved for. This procedure was continued until curves 
of l/p and (R.P.D)1/2 plotted against g passed through g = 0. The 
value for k was then cheinged and the procediire repeated. After suffi- 
cient variation of l/p and k, curves could be established of l/p 
against 2k-j^ = 2k(R.P.n)l/2 for particular values of the other param- 
eters M, g, f, and Pq/p* 


RESULTS AND DISCUSSION 


In the preceding sections a method of flutter analysis has been 
developed for a two-dimensional flat panel or plate held at its leading 
and trailing edges. The variables in the analysis are the number of 
modes or degrees of freedom the panel is assumed to have, Mach n\miber 
(greater than 1.0), l/p, 2ki = 2k(R.P.D)l/2, g, f, po/p, and edge 
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fixity. The analysis conveniently yields stability bovindaries in terras 
of l/^l and which are vised as the coordinates of most of the fig- 

ures presented. These two parameters are given in terms of the properties 
of the panel and svgiersonic stream by 


(30) 


where a is panel density, q is dynamic pressure, and is the 

first-mode eigenvalue given in table I for clamped and pinned edges. 
Inasmuch as the various parameters in the analysis contain implicitly 
the panel properties (E, a, v, and t/c), axial force F., air density, 
and speed of sound, the effects of varying these implicit properties can 
be obtained only by cross -plotting. 

In the present paper, some effects of the number of modes used in 
the analysis are studied by using two and then four modes of the clamped- 
edge panel with selected values of M, g, f, and Pq/p* addition, M 

g, f, and Pq/p varied in order to study their effects. To a lesser 

extent the pinned-edge panel is investigated for comparison with certain 
clamped-edge results . 

The following table lists the conditions for which stability bound- 
aries are given: 



Mach 

number , M 

Degrees of 
freedom 

Structural damping 
coefficient, g 

Tension 
parameter, f 

Density ratio, 
Po/P 

Clamped edge 

1.5 

2 

0 

0 

0 


1.0 

h 

0, o.cx)5, 0.025, o.o5> 0.05 

0 

* 0 

0 

0, 0.1, 0.5, 1-0 

0 


2 

0, 0.002, 0.00375, 0.05 

0 

0 

1.56 

2 

0 

0 

0 

k 

0 

0 

0 

Pinned edge 

CM 

0, 0.005, 0.00*175, 0.05 

0 

0 


The results are first grouped according to Mach mmiber and are later 
simimarized and compared. 
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Results for Mach Niunber of 1.3 

Effects of two and four modes .- Figirre 2 gives the results for the 
clamped- edge panel for the simple case of two degrees of freedom (first 
and second modes) with g = f = PqJp - 0* The abscissa is the stiffness 

0 >i c . 

parameter — — = 2k and the ordinate is the mass ratio 1/p. An ordinate 

of zero represents the limiting case of p = 0, or, in other words, the 
plate is vibrating in a vacuum. The two solid curves are the first- and 
second-mode stability boundaries as indicated. It was established by 
application of the Nyquist criterion (see, for example, ref. l4) as well 
as by interpretation of structural-damping resxats that the region to the 
right of the first-mode bomdary is stable, whereas the region to the 
left is unstable; furthermore, the region within the second -mode boundary 
is doubly unstable as indicated (lanstable with regard to both boundaries) . 
Values of the reduced frequency k are indicated along both curves. The 
points at which the curves cross the abscissa correspond to vibration in 
a pure normal mode (flutter at the limiting condition of p = O). 

It can be seen from equations (30), by taking the product of l/p 
and 2ki, that a specified panel material, air density, and speed of 
sound are represented by a hyperbola such as the dashed curve of figure 2 
with the panel thickness-chord ratio t/c increasing to the right. The 
intersection of the hyperbola with the stability boundary fixes the value 
of t/c for neutral stability. Thicker panels are stable and thinner 
panels are unstable. (The particular hyperbola shown is for aluminum 
panels in air with standard sea- level properties. For denser panels or 
less dense air, the hyperbola would be below the one shown.) 

Some effects of the number of modes in the analysis were studied by 
including the first four normal modes, and the results are shown as the 
solid curves of figure 3. The dashed cvirves are the results for two modes 
from figure 2. With the addition of the third and foxorth modes, the 
first-mode boundary is moved very slightly to the left (except where it 
crosses the abscissa) and is still the critical or decisive stability 
boundary separating the stable from the xmstable region. The second- 
mode boundary is also only moderately affected. Within the already 
unstable region there now exist third-mode and fourth-mode boundaries 
which are closely analogous in appearance and character to the first- 
and second-mode boundaries, respectively. The unstable region is divided 
by three of the boundaries into regions of different degrees of insta- 
bility as indicated by the numbers in parentheses ranging from (l) to (4). 
(The points at which the various boundaries cross the abscissa have the 
same significance as before.) These results indicate that two modes give 
a decisive stability boundary which is a close approximation to that for 
a large nvimber of modes, at least for the conditions g = f = pQ^p = 0 

and M = 1.3* 
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Effects of structural -damping coefficient g «- Flgiires 4(a) and 4(b) 
show the first-mode and second-mode boxmdaries (from a fovir-mode analysis) 
for various values of g (taken to be the same for all modes) . Third- 
and fourth-mode boundaries are affected by g in a manner similar to 
that of the first and second modes, respectively, and are not shown. 

The second-mode boundary of figvire 4(b) vanishes completely when g 
becomes slightly greater than O.O25, and for all positive values of g 
it remains in the vmstable region to the left of the first-mode boimdary. 
Included in figure 4 is the dashed hyperbola from figure 2. Since the 
thickness -chord ratio t/c decreases to the left in the figure, the 
abscissas of the intersections of the hyperbola with the stability bound- 
eiries in figure 4(a) show the proportional reduction in thickness reqviired 
to prevent flutter as g increases. 

Effects of tension .- Tension has a marked effect on the stability 
boundaries, as shown in figures and 5(t>) • Figure 5(a) shows the 

pertinent segments of the first-, second-, third-, and fourth-mode sta- 
bility boundaries for g = 0 eind for the three values 0.1, 0.5> and 1.0 
of the tension parameter f. As f Increases, all the boundaries move 
to the left, and the thickness required to prevent flutter is decreased. 
Futhermore, as f increases, the first-mode bovindary moves to the left 
more rapidly than the higher mode boundaries so that the rightmost bound- 
ary, separating stable from unstable regions, is one of the higher mode 
bovmdaries. For example, for f = 1.0 in figure 5(a), the third-mode 
. boxmdary is farthest to the right. This trend is not surprising since 
application of tension to the clamped-edge plate causes the largest per- 
centage increase in the first natural frequency, the next largest in the 
second natural frequency, and so on. Thvis, it appears that the Inclusion 
of only two modes in a flutter analysis may not be sufficient when the 
plate is subjected to tension. Inasmuch as the stiffness parameter 2k^ 

and the tension parameter f are both based on the first natural frequency 
of the panel without tension, the shift of the boundaries is due solely to 
the tensile force F. 

Figure 5(ti) shows segments of the first-, second-, third-, and 
fo\irth-mode boundaries with f = Ovl, 0.^, and 1.0 for g = 0.005* By 
comparing figures 5(^) and 5(a) it can be seen that g has a marked 
effect for the smaller values of f but its effect diminishes as f 
increases . 

Effects of still air below panel .- The one remaining parameter to 
be considered at M = I.5 is Pq/P^ ratio of the density of the still 

air below the panel to the density of the supersonic stream above. In 
the preceding results this ratio was zero. The effect of increasing 
PqJq to 1.0 will now be examined. For the sake of simplicity and con- 
venience, only a two -mode analysis is made. Effects of structural damping 
and tension are also included. 
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Figure 6(a) shows first- and second-mode houndaries for Pq/p = 1 
as solid curves and, for comparison, the dashed boundaries for Pq/p = 0 

from figure 2. Just as with the other results, the points where the 
boundaries cross the abscissa correspond to pvtre-mode resonance in a 
vac\xum. At these crossings the Imaginary part of C nti passes through 

zero. This imaginary part is a measure of aerodynamic damping. In the 
previous calculations Cnn consisted only of Inn^ whereas, for 
p^yip = 1.0, Cnn also contains Inn + ^nn (see eq. (23)). By compar- 
ison of the solid and dashed cvtrves on figure 6(a) it can be seen that, 
as a consequence, the first-mode boundary has moved to the left by about 
20 percent but the second-mode boundary has changed relatively little. 


Such an effect of still air might be expected since, for the same 
maximum panel amplitude, a first-mode vibration radiates into the still- 
air region a greater amount of energy per cycle than does a second-mode 
vibration. (With regard to the radiation of sound from a piston in a 
plane wall, specifically for the case of a piston with nonrigid face, 
p. 336 of ref. 15 gives the result that, at frequencies which are small 
compared with the ratio of the speed of sound to 2 jc times the piston 
radius, the pressure on the piston is approximately uniform and nearly 
proportional to the average velocity of the piston. Since the average 
velocity of the second mode and all other antisymmetric modes is zero, 
the pressure due to these modes is nearly zero and, accordingly, almost 
no work is being done on the still air.) 

As can be observed in figure 6(a), the first-mode boundary has moved 
to the left of the second-mode boundary in the region of small mass ratio; 
in this region the second-mode boundary becomes critical. 

Figure 6(b) shows the effects of structural damping on the first- 
mode boundary, which for g = 0 is shown more completely in figure 6(a). 
Cvtrves are included for g = 0, O.OO5, 0.05, and O.O5. For values of g 
larger than about 0.025 the second-mode boundary vanishes as it did pre- 
viously with Pq/P = 0 in figure 4(b) , and only the first-mode boundary 

remains. The dashed hyperbola for alvimlnvim at sea level is included in 
figure 6(b) , and it can be seen, that a plate with zero structural damping 
would have to be about 30 percent thicker than one with g = O.05 in 
order to prevent flutter. 

Figure 7 shows the effects of the tension parameter f for pp/p = 1 

and g = 0. Both first- and second-mode boundaries are shown for f = 0, 
0.1, 0.5, and 1.0. In this case, just as with pp/p = 0, tension causes 

a marked reduction in the thickness required to prevent flutter. Futher- 
more, if more than two modes had been included, tension would have resxilted 
in a higher mode bovindary farther to the right than the cxirves shown for 
the higher values of f. 
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Resvilts for Mach Number of \f 2 

Clamped-edge panels .- Figirre 8(a) shows the stability boundaries 
obtained from a two-mode analysis for clamped-edge panels at M = \[2 
with f = pq/p = 0 for various values of g. Included in the fig\ire is 

the dashed hyperbola appropriate to this Mach number for aluminum panels 
in sea-level air. From a comparison of figure 2 and the curves of fig- 
ure 8(a) for g = 0, the first-mode bdmdary of figure 2 appears to have 
moved into the positive mass-ratio region and the second-mode boundary 
appears to be moving toward the negative mass-ratio region. Such is the 
case, but, inasmuch as the flutter frequencies on the upper boimdary of 
figure 8(a) are about midway between the first and second natural fre- 
quencies, this bomdary can now be referred to only loosely as a "first- 
mode" boundary. The lower boxmdary is still readily identified as a 
second-mode boundary and the intersection with the abscissa corresponds 
to vibration in a pvure second mode. 

In contrast to the situation at M = 1.5> the second-mode boundary 
for g = 0 is now decisive for panels represented by the dashed hyperbola. 
Values of thickness to the right of the second-mode boundary are stable 
and, in addition, a small range of thickness values is stable between 
the upper and second-mode boundaries. 

The c\irves in figure 8(a) for positive values of g show that the 
region of instability within the second-mode boimdary is reduced for 
small values of g (as for M = I.5) and vanishes when g is slightly 
greater than O.CX)575^ tut that small values of g increase the region 
of instability associated with the upper boundary. This effect of g 
on the upper boundary is in marked contrast to its effect on the first- 
mode boundary at M = l.J. (See fig. 4 (a).) The differing effects of 
structural damping at M = 1.5^ M = \/ 2 , and M = I.56 are considered 
further in the section on "Variations With Mach Number." 

Pinned-edge panels .- In order to indicate effects of edge fixity, 
bovmdaries are shown in figure 8(b) for conditions identical to those 
of figure 8(a) except that the edges are pinned rather than clamped. 

The boundaries for g = 0, which are given incompletely in reference 5 
and thereby lead to the conclusion that only a small range of panel 
thickness is stable at M = \l~ 2 , have been extended to higher frequencies 
with the result that siifficiently thick panels are also fovind to be stable. 
The effect of structural damping on both boundaries in figure 8(b) is very 
similar to that in figure 8 (a). The dashed hyperbola appropriate to 
pinned-edge aluminum panels in sea-level air is included in the figure. 

The hyperbolas of figures 8(a) and 8(b) are located differently because 
of the different values for the first-mode eigenvalue K2^ of equation (50) 

for pinned and clamped edges. (See table I.) From the intersection of 
the dashed hyperbolas with the stability boundaries in the two figxores, 
it can be determined that a pinned-edge panel must be somewhat thicker 
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than a clamped-edge panel in order to be flutter free but not nearly as 
thick as might be expected from a simple comparison of the first natural 
frequencies . Values of the reduced frequency k are indicated along each 
of the boundaries of figure 8. 

Based on what occurred at M = l.J (see fig. 5) > there is the pos- 
sibility that for g = 0 the fourth-mode boundary from a four -mode analy- 
sis would alter the stability picture in both figures 8(a) and 8(b) in 
the relatively unimporteint narrow range of stability between the bound- 
aries shown for g = 0. This minor effect of the fourth-mode boundary 
is expected to disappear for values of g greater than about O.OO 5 and, 
therefore, a four-mode analysis was not made for this Mach number. 


Results for Mach Number of I .56 

Effects of two and fovir modes .- As in the case of M = 1.5, stability 
boundaries were obtained first for two and then for four degrees of free- 
dom with g = f = Pq/P = 0* These boundaries appear in figure 9 as dashed 
curves for two modes and solid curves for fovir modes. Values of the 
reduced frequency k are indicated along the boundaries. The stable 
region is again to the right, and on the left the degree of instability 
is indicated in parentheses for the four -mode analysis. 

The two-mode results in figure 9 continue the trend noted in the 
preceding section from comparison of the curves of figure 2 and those of 
figure 8(a) for g = 0. The second-mode boundary has moved entirely into 
the negative mass -ratio region. The upper boxmdary has moved higher in 
the positive mass-ratio region, and the flutter frequencies along it, 
which in figure 8(a) were midway between the first and second natural 
frequencies, are now closer to the second. For this reason the upper 
boundary, which was loosely identified as a "first-mode" boundary in the 
discussion of figure 8(a), will now be referred to as a "second-mode" 
boundsiry. 

A further point of difference between the results at M = 1.3 and 
M = 1.56 is that the addition of the third and foixrth modes at M = I.56 
shifts the decisive stability boundary to the left by about 10 percent, 
whereas at M = 1.5 the shift is insignificant. (Compare figs. 5 and 9*) 
Although this shift indicates that the two-mode result is not well con- 
verged, the two-mode boundary is conservative; that is, it requires a 
greater thickness to prevent flutter. (As with two modes, when four modes 
are used, half of the stability boxindaries fall in the negative mass- 
ratio region. ) 

Effects of structural damping coefficient g .- No curves are shown 
to indicate effects of structural damping at a Mach number of I. 56 , the 
reason being that, for moderate values of the coefficient g, ranging at 
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least up to 0.05, the stability boundaries fall virtvially on top of those 
for g = 0. The major effect of structural damping is a moderate change 
in flutter frequency. 

Effects of tension and of still air belov panel .- Effects of tension 
have not been determined, but tension is expected to have essentially the 
same favorable stiffening effect at all Mach numbers as at M = l.J. The 
effect of still air behind the panel has also not been determined, but 
this effect is expected to be less than at M = I.3 for two reasons: 
First, the air beneath the panel acts primarily as ein energy absorber and 
one means of energy absorption, structural damping, has been found inef- 
fective in shifting the stability bovindaries. Second, on the decisive 
boundary the flutter mode appears to be predominantly the second natural 
panel mode, and it was fovind tliat at M = 1.3 the second-mode boundary 
is changed only slightly by increasing pq/p from 0 to 1. 


Variations With Mach Number 

The foregoing res\xlts have been presented for particular Mach num- 
bers. In an effort to clarify some of the anomalies that have been noted 
in these results, figvires 10 to 12 are presented. Figure 10 shows the 
panel thiclmess -chord ratio required to prevent flutter as a function 
of M for clamped-edge panels with g = f = p^^p = 0. The curves apply 

to aluminum panels in standard sea-level air. The values at M = 1.3> 
f2, and I.56 were obtained from figures 2 and 8(a) and the two-mode 
results of figure 9- The shape of the curves between these known points 
is estimated. The stable region is above or to the right of the shaded 
boundaries . 

The boundar-y which is labeled "first-mode" on one end and "second- 
mode" on the other has flutter frequencies which progress from slightly 
above the first natural frequency to somewhat below the second natural 
frequency as the Mach ntimber is increased. (See previotis discussions 
concerning figs. 2, 8(a), and 9-) The boundary labeled "second-mode" 
has flutter frequencies slightly below the second natiiral frequency 
throughout . 

Figvure 10 shows the second-mode stability bovmdary to be decisive 
in the Mach number range from slightly above I.30 to slightly above )[2. 
As the structural damping g is increased from zero^^ the second-mode 
boundary shrinks to the left leaving the "first-mode" — "second-mode" 
bovmdary decisive throughout the range of M shown. For example, for 
a value of g slightly greater than O.OO38 the second-mode boundary does 
not exist at U = \[2 (see fig. 8(a)), and for a value of g slightly 
greater than 0.025 it does not exist at M = 1.3 (see fig. i<-(b)). 
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These effects of g on the second-mode boundary are illustrated In 
figure 11, which contains cross plots of g against t/c obtained from 
the intersections of the dashed hyperbolas (for aluminum panels in sea- 
level air) with the boundaries for constant g such as shown in figures k, 
8(a), and 9- Figure 11 also shows that an increase in g from zero wotild 
ca\ise the "first-mode" — "second-mode" boundary of fig\are 10 to drop mark- 
edly at M = 1 . 3 , rise slightly at M = \/2, and remain essentially 
unchanged at M = I.56. The ratios of flutter frequency to the first nat- 
ural frequency (o/(Oi are indicated for each of the crossings end tend to 


show more clearly the connection between figures 10 and 11. (Values 
of co/o)! near which is approx i mately 2.76, are associated with 

the second-mode boundary and values between 1.0 and 2.05 associated 
with the "first-mode"— -"second-mode" boTindary. ) A complete xuiderstanding 
of the manner in which the curves change character and position with Mach 
number, particularly between M = 1.5 and M = \/2, reqiilres more cal- 
culation than presented herein. 


Figure 12, which has the same ordinates as figure 10, is presented 
for the ptirpose of summarizing some effects of all the parameters investi- 
gated in the present paper. The results shown are based on two modes, 
except in the case of tension where only four -mode results are known. 

The results again apply to almimma panels in sea-level air. The fig- 
ure shows as a solid curve the shaded bovindary from figure 10 for clamped- 
edge panels and as a short-dash curve the effect on this boundary of 
increasing g from 0 to O.O5. The third (long -dash) cijrve is for plnned- 
edge panels with g = 0, the value at M = 2 having been obtained from 
reference 2. The points at M = /2 were obtained from figure 8(b) and 
the upper (second-mode) ctirve was patterned after that for clamped-edge 
panels. As a matter of interest, points are included in figure 12 at 
M = 1.5 foi” clamped— edge panels and indicate the effects of tension 
(f = 0.5) and of still air below the panel (pq/P = 1-0) for 8=0 

and g = 0.05. 


Some effects of the various parameters can be assessed from figure 12. 
The overall resiilt is that t/c is highest in the low supersonic Mach 
number range and suggests that this range is the more critical from a 
design standpoint. Structural damping is seen to have a large favorable 
effect near and below M = V2. Although rather influential at M = 1.5^ 
the still air below the panel is expected to have less effect at M = \/2 
and 1.56. Tension, which is seen to have a large favorable effect at 
M = 1.5, is expected to be similarly effective for all Mach numbers. In 
this connection, it might be mentioned that one means of producing tension 
is by a static-pressure difference between the upper and lower panel sur- 
faces, particularly for the case where the panel leading and trailing edges 
are prevented from moving toward each other. A comparison of the results 
for the edge fixities, pinned and clamped, is of interest because the 
edge fixity of actvial panels falls somewhere between. 
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Comparison With Other Work 

In reference 2 the conclusion is reached that all panels, regardless 
of thickness, are unstable for supersonic Mach nvntibers less than \/2. This 
result eind the more plausible results of reference 2 for M > V2 are based 
on air forces expanded to the first power of the frequency of oscillation. 

In reference 5 the necessity of including higher order frequency terms for 
Mach numbers near \J2 is pointed out, and stability boundaries, based on 
exact linearized unsteady air forces, are presented for M = \/2 and M = 2. 
One boundary is obtained at M = 2 which agrees well with the comparable 
resxilt from reference 2, whereas two boundaries are obtained at M = V2. 

The boundaries at M = \/2, because they are not carried high enovigh in 
frequency, are interpreted in reference 5 as showing that stability is 
possible only for a small range of peinel thickness at this Mach number. 

The resvilts of reference 2 for K < \f2 are not questioned in reference 5* 

In the present paper, stability boimdaries are computed for M = \/2 
and for Mach manbers above and below this value (namely, M = 1.3 and 
M = 1.56). In view of the findings of references 2 and perhaps the 
most noteworthy result of the present investigation is that, for the Mach 
numbers treated and probably throughout the supersonic range, sufficiently 
thick panels are stable. 

In references 2 and 5 and the present paper, M = \/2 appears as a 
transitional value. The transition is evidenced herein by the contrasting 
behavior of the stability boundaries at M = 1.3 and M = I.56. Some 
understanding of why a Mach number of \^2 is transitional can be had by 

examining matrix eqijation (26). The term — ~ |zn| equation (26), 

Pkp3 *• J 

being the entire first-order frequency contribution to the damping, is 
dominant at low frequencies. This term appears to control the slope, at 
low frequencies, of the eventually decisive stability boundary and changes 
sign as M passes through \/2. When M < \/2, the sloi)e is negative for 
low frequencies, but as the frequency increases the slope eventually 
becomes positive because of the higher order frequency effects (for example, 
in fig. 2) . Because only first-order frequency effects are included, in 
essence only the beginning portions of the stability boundaries for 
M < \/2 are obtained in reference 2, and, as a consequence, the conclusion 
is reached that eill panels are unstable below this Mach nvimber. For 
M > \/2 the slope of the decisive stability boxmdary starts out positive 
and becomes more so as the frequency increases. (See fig. 9*) If aspect 
ratio were included in the present treatment (by considering three- 
dimensional rather than two-dimensional panels) , a reduction in aspect 
ratio would probably tend to eliminate the initial negative slope of the 
eventually decisive stability boundary for M < \f2 and increase the 
initial positive slope for M > \/2. This effect of aspect ratio is 
expected becavise, in general, a reduction in aspect ratio results in an 
increase in aerodynamic damping with a consequent enleirgement of regions 
of stability. 
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In reference 7, which treats a different problem (namely, an infinite 
two-dimensional panel on equally spaced supports) , the result was also 
obtained that somewhere in the s\ipersonic Mach number range a panel will 
flutter regardless of its thickness. The conclusion was reached that 
stability is not possible at supersonic Mach nvmibers less than about 1.25 
and that at higher Mach numbers a sufficient increase in thickness will 
always render a stable panel unstable. However, it was observed that 
over a large portion of the predicted region of instability the flutter 
was of an extremely mild character, since a large number of oscillations 
were required to double the amplitude. With the hope of eliminating the 
large region of mild instability, small amounts of viscous damping were 
included. Contrary to expectations, thick panels remained unstable for 
the example given at M = 1.8. 

As part of the viscous -damping investigation, the results were inter- 
preted so as to determine regions of stability and instability. As shown 
in figure 10 of reference 7, an apparent conflict with the results of 
Nyquist diagrams was found. (The Nyquist diagram concept is used in 
general in reference 7 for determining stability. ) This conflict is based 
on the assumption that most investigators interpret structural-damping 
resvilts according to the concept that removal of damping tends to destabi- 
lize. This assumption is incorrect, however, and no such simple criterion 
holds true for interpreting struct\aral-damping results. A feat\are to be 
noted in the example chosen in reference 7 to illustrate the apparent 
conflict is the existence of infinite singularities in the air forces at 
the end points of the boundaries (g = 0, 0.01, and 0.05) on the right 

in figlire 9 of the reference (designated type B loci therein). By way 
of explanation, such singularities occur in the linearized-flow treatment 
becaiise a traveling wave of panel deflection is moving at a speed cor- 
responding to M = 1 relative to the air above or below the panel. 

In preparing the present paper, the question of stability was inves- 
tigated by means of both the structural-damping concept and the Nyqviist 
diagram concept. The structural-damping resxilts in every case agreed with 
the Nyquist diagram results. 

Incidentally, in using the Nyquist concept, knowledge of the aero- 
dynamic forces- for all frequencies from minus infinity to plus infinity 
is required.' Thus, the concept is not applicable, in general, when the 
air forces are approximated by a few terms of a power-series expansion 
in the frequency of oscillation. 


CONCLUDING RByiABKS 


A Rayleigh type analysis involving chosen modes of the panel as 
degrees of freedom has been used to treat the flutter of a two-dimensional 
flat panel supported at its leading and trailing edges and subjected to 
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a middle -plane tensile force. The panel had a supersonic stream passing 
over its upper stirface and still air below. The aerodynamic forces due 
to the supersonic stream were obtained from the theory for linearized 
two-dimensional unsteady flow and the forces due to the still air were 
obtained from acoiastical theory. The still air beneath the panel was 
treated on the assumption that the still-air reservoir extended to 
infinity. Accordingly, once acoustic energy was radiated into this region, 
none of it was ever reflected. Such a situation is, of course, not the 
same as for a panel on a closed body but represents a first approxima- 
tion for many practical cases. 

In order to study the effect of increasing the nimiber of modes in 
the analysis, two and then foiu: modes were employed. The modes -used were 
the first four nat\oral modes of the panel in a vacuvm with no tensile 
force acting. The analysis included the variables: Mach number M, struc- 

tural damping, tensile force, density of the still air, and edge fixity 
(clamped and pinned). For certain combinations of these variables, sta- 
bility bomdaries were obtained which can be used to determine the panel 
thickness required to prevent flutter for any panel material and altitude. 

In contrast to some previous panel flutter investigations, the results 
of the present paper show that sufficiently thick panels are flutter free 
for the Mach numbers treated and suggest that this is true throughout the 
supersonic speed range. The low supersonic Mach numbers were found to be 
most critical from a design standpoint in the range examined (from M = 1.5 
to M = 2 , 0 ). Tension was shown at M = 1.5 to have a marked favorable 
effect (also expected at all Mach mmibers) in reducing the thickness 
required to prevent flutter, and it was pointed out that one means of 
producing tension is by a static pressxxre difference between the upper 
an «^ lower sirrfaces of the panel. Small amovints of structural damping 
were found to have a pronounced beneficial effect near and below M = \/2 
and essentially no effect at M = I. 56 . In the neighborhood of M = \/2 
a small change in either Mach number or structxiral damping was foimd to 
cause an abrupt change in the thickness required to prevent flutter. At 
M = ^2 a pinned-edge panel must be somewhat thicker than a clamped-edge 
panel in order to be flutter free. Still air below the panel was taken 
into account only at M = 1.5 and was shown to have a moderate beneficial 
effect. For M > V2 the still air is expected to have little effect 
because for this Mach number range the flutter mode is predominantly the 
second natural mode, which radiates very little energy into the still air. 
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It may be worthwhile to mention that two factors not investigated 
in .the present paper, aspect ratio and panel buckling, will probably 
have strong effects on the thickness required to prevent flutter. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 20, 1955 • 
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APPENDIX A 

SOLUTION BY MEANS OF LAPLACE TRANSFORMS 


In a recent paper (ref. 6) , a procedure is outlined for obtaining 
by means of Laplace transforms the exact solution for the flutter of a 
two-dimensional membrane which is subjected to a supersonic stream on 
one side and stagnant air on the other. This solution is called exact. 
Inasmuch as the equation of motion for the system is solved directly 
without any limitation being imposed on the mode shape or frequency of 
flutter. Reference 6 also mentions that pure bending of a plate and the 
more general case in which plate bending and membrane stretching are com- 
bined could be treated in the same manner. The present paper treats the 
latter case; namely, the flutter of a panel (plate) acted on by a middle- 
plane or axial force, such as tension or compression, less than the buckling 
load. In the body of the paper this problem is solved by the generalized- 
coordinate approach, and the coordinates used are the normal modes of the 
panel with no axial force acting. In this appendix, for the sake of com- 
pleteness, the solution to the same problem is derived by means of Laplace 
transforms to the point where numerical calculations can be made. The 
feasibility of applying the Laplace transform solution is examined, but 
no nvimerical results are obtained. 

The integrodifferential equation to be solved is given by equation (2) 
which, upon substitution of the expression for p(x,t) obtained from 
equations ( 5 ) to (9)^ may be written as 


^ Z"" - o^m^Z - J1 Z" + 

c2 c 


r wjoijx- 


dx J 


0 




1 ^ 

2 P 


i2k f wj|)lu(x-0d| + 
^0 

2k r Wi(|)l^ 

^0 


(x-|)d| 


(Al) 


or alternatively as 


^ z"" - - Z- Z" + Ji 

c^ c2 c 


Wu(0)lu(x) + 




1 1 ] 

i2k / Wj(0lz(x-0d| > 

0 JJ 


= 0 


(A2) 



52 


NACA TN 3^6^ 


where 


wjx) = Z'(x) + i2k Z(x) 
w^(x) = i2k Z(x) 


Ijx) = e-i^Jo ^ X 


( 2 ), 


Ij(x) = ''(2 M|x|) 

Dividing equation (A2) by mcu? yields 


aZ"" - &Z" - Z + - - 


2 P 
where 


i2k jT wj(|)lj(x-|)d| 


+ i2k\wu(|)d^ 


Wu(0)lu(^) + / ■^UV'^'S^ITT -^“-|"uv 

J 0 


> = 0 


(A5) 


= J:- a 
h\^ I 4 

> V 


6 = ^ . - m =m 

2 2\(i^ 

m^c o^\ / 


The quantity coq in the formulas for a and 6 is the first natural 

frequency of the plate vibrating in a vacuum with no axial force F 
acting and Kq is the associated eigenvalue. (See table I.) In the 

case of the membrane, ouq would be the first natural frequency of the 
membrane, a would be zero (D is negligible for the membrane), and f 


would be I - 

lit 
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Applying the Laplace transform 


L {z(x)J = Z(s) = J e"®*Z(x)cbc 


to equation (A3) (using in the process pairs 4 and 'J, p. 294, ref. 16) 
yields the transformed problem in the form 


a^s^Z - s^Zq - s^z-^ - szo - Zx) - bfs^Z - sz, 


.2 - - 5 ^ 


^s^Z - szq - Zj^j - 


Z + 


— f— ^s^ + i4ks - 4k^ Z - (s + i2k)zQ 1 ^( 3 ) 

z u \P — — 


4k^i 


1 ^ 4k2L 

2 P 


f Z(|)lj.(x-|)d| 
'■>0 


= 0 


(a4) 


where Zq = Z(0), Zj^ = Z'(0), Z 2 = Z"(0), z^ = Z"'(0) , and Iu(s) is 
the Laplace transform of Iu(x) • The Laplace transform in eqviation (a 4) 
involving Pq^P as a multiplier is the contribution of the perturbation 

pressure on the lower surface of the panel. Unfortimately, this transform 
does not appear to be .obtainable here where the deflection Z is unknown. 
In the body of the paper the effect of including the air below the panel 
is foxmd to be moderate at a Mach number of I .3 and reasons are given why 
this effect is expected to be even , smaller at the higher Mach numbers 
investigated. In view of these facts and in view of the difficulty of 
handling the lower-surface term in equation (a 4), this term will be 
omitted in the rest of this appendix - that is, treated as if pq were 

zero . 


Equation (a4) 


can therefore be reduced to 


^os4 - 6 s 2 - 1^ + £(s + i2k)2i^(s) Z(s) = 


s^z, + sz^ + z 


>)- 


6 Zn 


(A5) 


where e 



and Zq has been dropped because it is zero for the 


boundary- value problems of the present paper. Thus the integrodiffer- 
ential eqxxation (A1) has been reduced to the algebraic equation (A5) • 
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Now by means of pair 11, p- 29^, and pair 55^ P* 298, of refer- 
ence 16 there is obtained 


iu(s) = 


(s 




-1/2 


(a6) 


Therefore, from equation (A5), after some algebraic manipulation 

^ M(s) N(s) - , 

Q(s) Q(s) 


(A7) 


where 


Q(s) 

M(s) 


= (■ 


os^ - 6s^ - 1^ 


= ^as^ - 6 s2 - 1^ 


(s . . (if 


(s + ioj)^ + 

\ M , 


- e2(s + i2k)^ 
a^s^Z]^ + SZ2 + - bzj^ 


N(s) = - €(s + i2k)‘ 


(s . io^)2 . I 


a^^s^Zj^ + SZ2 + - 


6zi 


In polynomial form the quantities Q, M, and N are 


10 


Q(s) = H q^slO-r 


(a8) 


r=0 


M(s) = 


H Z + Z2 £ nv(2)s'?-^ t .3 51 nv(3)s«-- (A9) 


r=0 


r^O 


r=0 


N(s) . ^ nr(l>s8-'' * ■£ nr(5)s'*-*' (AlO) 

r=0 r=0 r=0 
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The coefficients of the various series are given in the following table: 


r 

Or 


«r( 2 ) 

mr(3) 

Br(l) 

n,( 2 ) 

ni.(3) 

0 

a? 

a 2 

a 2 

a 2 


-ae 

- 0 € 

1 

12 cu »2 

1250,2 

125i2 

iaai 2 

-i25«(L.g) 



2 



^ 

.^05 + 02 ^) 




5 


.11(5x8 

-125o8 

-laixs 

lacfa/i + Ei) ♦ a, el:^] 
[ \ M^/ W* _ 


i5«8!;^ 

y* 

k 

S 2 . aa + oS 2 !^ 
((2 

52.a+aaeM 

^ + 056 ^ 

o5 SM 


m 6 


5 

125(62 - 2a) 

125(82 - a) 

-125a 

-12uxx. i 


0 


I 

25 - €2 - (82 . 2a.)^ 

5 . (82 - a)l!^ 


m2 

8e£^ 



1 

i(l*5B - at2c2) 


0 i 





1 

1 . 28 + 24k2c2 

m2 

IQ 




’ 


i 

1(25 + 32k^€2) 







K) 




1 

i 

i 





The exact inverse transform of equation (A 7 ) requires the determina 
tion of the roots of Q(s) (eq, (A8)). Since Q(s) is a tenth-order 
polynomial, its roots can be solved only approximately for specific 
values of the coefficients q^. An alternative procedure is to expand 

the quantity [q(s)] in a Maclaurin's series (a procedure used in 
ref. 17) , with the result that it may be expressed in the form 


where 


1 

Q(s) 


T ^ T 

- 

slO ^ sn 


qoTo = 1 


(All) 


10 

Qo^n “ " — ^r^n-r (^ = 1) 

r=l 

and T with a negative subscript is to be interpreted as zero. 
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When the series expansion for _Cq(s)J” (eq. (All)) is substituted 
into equation (A7), the transform Z(s) becomes the sum of infinite series 
with terms of the two distinct types 

A 


and 




where m is a positive integer. The inverse Laplace transform of the 
first type of term is (see pair 5, p. 295 of ref. 16) 



AxDI-1 

(m - 1): 


and of the second is (see pair "J, p. 294 of ref. 16) 


(A12) 


L 


-1 




(x - |)”'"llu(0<i| 


(AI5) 


where Iu(x) is defined following equation (A2) . 


Upon substituting equations (A9)^ (AIO), and (All) into equation (A7) 
and using equations (A12) and (AI3) to obtain the inverse transform of 
the resultant expression, Z(x) is given by 


Z(x) = h3,(x)z2^ + h2(x)z2 + h^(x)z^ 


(A14) 


where 


” 8 _ T „ ( 1 ) vn+r+l « ^ m „ ( 1 ) px 

n =0 r =0 (n + r + l)i n =0 r =0 (n + r + 5 ).' '-^0 

" .J— m m ( 2 )_n+r +2 " 5 m n (^) r>x , 

= Z Z 7 “ -, -r r * H J1 r ° M ' f - E) Iu(E)dE 

n=0 r=0 (n + r + 2 ) . j^q p-Q (n + r + 4), Jq 

^ /-.x 


- 8 m (5)vn+r+3 

= Z_ ^ 7 

n=0 1^0 (n + r + 5 ) . 


Tnnr^5) 




t\ 


n+r+*5 
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In deriving equation (AIU) only one boundary condition - namely, 

ZQ = Z( 0 ) =0 - has been used thus far. In order to obtain the solu- 
tion for a plate restrained in a particular manner, it is necessary to 
impose three additional boundary conditions. These additional conditions 
for the plate with pinned and clamped edges are given in equations ( 5 ) 
and ( 4 ), respectively. By their use, one of the terms of equation (AlU) 
is eliminated and two homogeneous equations in the two remaining unknown 
Zj^'s are obtained. The borderline condition of harmonic oscillation, 
or the point at which flutter occirrs, is obtained by setting the deter- 
minant of the coefficients of these equations equal to zero. Thus, the 
flutter determinant for the pinned-edge plate is 



hi(l) 

^5(1) 


hl"(l) 

f— 1 

and for the clamped-edge plate is 



h2(l) 

h^(l) 


hp’d) 

hj ' ( 1) 

where the determinant elements are 

given by 


(AI5) 


(ai6) 


° S S (n + r + 1 ); S I 5 (n + r + 3 ); Iq lu( 0 <i| 

hl”(l) = 1 Z -^ 2 :^ ^ Z Z ' ■ . f (1 - 

n?tl (n - 1 ), n =0 r=l (n + r - 1 )! n =0 r =0 (n + r + l ) ^0 

Ml) = z z , * z f: (1 - 

n =0 1^0 (n + r + 2). n=0 r=0 (n + r + 4 ) : J q 

°> oo 5 qin( 2 ) pi 

h2-(i) = Z Z - Z Z / (1 - 

n =0 r =0 (n + r + 1 ). n =0 i^O (n + r + 5)1 '^o 


hjU) = 


^ Tnny(3) ^ ^ ^ Tl 

rt=0 r=0 (n + r + 3) I i^O r=0 (n + r + 5 ) I "^0 


-/ (1 - 
1 J n 


h 3 '(i) = Z Z -^2^ . Z Z /^l - 0 --\( 0 d. ' 

n=0 r=0 (n + r + 2) . n=0 i^O (n + r + 4) J '^0 


h 3 "(l) = Z II , ^ Z 7 -- — / (1 - 

It=0 1^0 (n + r + 1) ; n=0 1^0 (n + r + 3) J ^ 0 


(5) 
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Each of the preceding elements contains integrals of the form 


Im(M,oj) 



l)“lu(l)d| 


(AI7) 


which can be written in terms of the Schwarz fmctions f^(M,ffi) (see 
ref. 9 or eq. (24)) as 


Im(M,aj) 


f: 


(- 1 )^! 

fx(M,oi) 

(m - a)ia: 


(A18) 


where 



Examination of the series in the elements of equations (AI5) and (A16) 
reveals that A of equation (AI8) ranges at least between 0 and n + 7 
and at most between 0 and n + 9- In order to obtain accuracy to four 
significant figures, at least the first eight terras of each series and 
the consequent ranging of A between 0 and 16 are probably required. 
Inasmuch as the Schwarz fvinctions f^ have been tabulated for only the 

first few values of A, the use of equation (AI8) would require the deter- 
mination of a rather extensive series of fA's- An alternative and per- 
haps more efficient procedure would be to evaluate directly the inte- 
grals as given in equation (AI7) rather than to resort to the expanded 

form in equation (a 18) . 

Attention will now be given to the solution of the determinantal 
equations (AI5) and (A16) . A method of solution for parameters that were 
sought in the generalized-coordinate approach of the body of the paper 

^that is, l/p and 2ki = 2k(R.P.a)^/^) will be outlined here. 

The elements of eqmtlons (AI5) and (A16) are complex functions of 
the five parameters M, k, 0, (with g = O) , f, and l/p.* The most 
difficult parts of these elements to evaluate are the integrals generi- 
cally represented by in equation (AI7) , which are functions of the 

parameters M and k. Therefore, a convenient method of solution would 
be to fix the parameters M and k and preferably f and vary the 
remaining parameters 'and l/p in the left-hand side (hereinafter 
referred to as A) of equation (AI5) oi" of equation (A16) . By varying Q 
and l/|i over sufficiently broad ranges, an indefinitely large niimber 
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of combinations of and l/|i which cause A to vanish could be 
found. Each combination would define a point on separate stability 
boundaries, such as those shown in figure 3* Each boundary could then 
be determined as completely as desired by varying k over a sufficient 
range and repeating for each chosen value of k the process of finding 
combinations of and l/p which cause A to vanish. 

As can be surmised, the numerical calculations would be extremely 
lengthy even apart from two other questions which arise; namely, which 
is the stable side of each boundary, and has the critical bo\indary been 
fo\ind which separates stable and unstable regions and thereby defines 
the thinnest panel that is stable? In the present paper, therefore, the 
stability boundaries shown in figvires 2 to 9 were calculated exclusively 
on the basis of the generalized-coordinate or modal approach. 
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APPENDIX B 

VELOCITY POTENTIAL AND RELATED INTEGRAL 
Velocity Potential 0 ^ 


The velocity potential given in equation (9), which applies to 

the lower sirrface of the two-dimensional panel shown in figure 1, will 
now be derived. The system consists of a panel of width c, which is 
part of an otherwise rigid surface of infinite extent, oscillating har- 
monically with stationary air extending to infinity below. Thus, over 

the panel the normal velocity on the lower surface is Wj = ico Z(x)e^’^, 
while over the rest of the plane w^ = 0. 

According to reference 10 the solution to this problem can be 
obtained from 



(Bl) 


where w^ is the given normal velocity at the element of area dS of 

the plane and 0 is the velocity potential at a point P which is at 
a distance r from dS. From equation (Bl) the velocity potential at 
the surface of the panel may be obtained, in terms of the coordinates of 
figure 1, as 





Z(|)d| 



-iiix-OW 

> = dy 

\/(x - 


(B 2 ) 


Upon making the substitution y = |x - 5 | cosh G, the integral with 
respect to y in equation (B 2 ) may be written in the form 


l(x-0 = 




dy = 


i^|x-5 cosh 0 


de 


(B 3 ) 


+ y 


-.00 


/(x - 1)^ 


~00 


e 
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By means of equation (11) on page 18O of reference I8, equation (B5) 
becomes 


I(x-I) = (b4) 

Substitution of equation (b 4) into equation (B2) yields 


If the .coordinates x and ^ are nondimensionalized by dividing by the 
panel chord c, the form for given in equation (9) is obtained. 


Integral 

The third term on the right-hand side of eqioation (22), namely. 


Pn(^) = ^- Zn(Ologe(kM|x - | | )d| (b6) 

P ^ Jo 


contains the singularity of the Hankel fimction in equation (20). 

As a first step in the evaluation of the improper integral in 
equation (b 6), let 


and 


I = i(l _ cos 0 

> 



cos t) 


\(t) 


00 



m=l 


Sm sin m^ 


(B7) 


(B8) 
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where 



2^U)s±n mi di 


In terms of equations (B7) and (b 8), eq\aation (b 6) 'becomes 


Pn(x) = ^ 

p 2it 


(B9) 


where 




m=l 




/ 

sin i sin m^ log„(- 

— Icos i - cos \lf| ) 

Jo ®VS 

! ' V 


(BIO) 


Taking the derivative of equation (BIO) with respect to f and making 
use of reference I 9 to evaluate the resulting improper integrals yields 


M. 

dijr 

= ^ S|jj jsin(m - l)\lr - sin(m + l)>|^ (Bll) 

^ m=l 


Integration of equation (Bll) gives 


^ fee Jo ^ - cos 


I(^) = I 


^|Si cos 2t + 


m=2 


Sm 


cos(m + l)t 
m + 1 


cos(m - l)\lf 


m - 1 


+ K (B12) 


The integration constant K in equation (B12) is determined by setting 
equal to jt/2 in equations (BIO) and (b 12) and eq\iating the two 
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resultant expressions. By so doing, it is found that 


K = log^ M 


By means of equations (B 12 ) and (BI3), equation (B 9 ) becomes 


P„(X) 


^0 1 A kM 1 „ cos(m + l)\lf 

p 4 |V ®4 2 71^^ m +1 


cos(m - l)\|f 
m - 1 


The mode shapes in equation (b 6) are now approximated by the 

finite sine series 


Zn(0 / Sr sin rx 
r=l 


The constants Sj^ in equation (B 14 ) are obtained from the expression 

following eq\iation (b 8 ), with the resiat that equation (B 14 ) can be 
written as . 


Pn(x) = Lj.(x) 

P ^ r=l 


Lj^Cx) = logg i COS 2 t 


Iy(x) = — cos(r + l)t i — cos(r - l)\|r (r > 2) 

^ r + 1 r-1 ' = / 
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and, as in equation (B7), 


X = i(l - cos \|f) 


The form given in eqiiation (B16) was used to obtain eqxiation (28). Of 
interest is the fact that only depends on k and M. The term 

is therefore comparatively simple to Include in equation (22). 

The coefficients s^. in equations (BI5) and (B16) for the first 

four modes of the plate with clamped edges are given in the following 
table : 



Mode 

1 

Mode 

2 

Mode 

3 

Mode 

4 

®1 

0.66615 

0 

-0.52917 

0 

82 

0 

0.64119 

0 

0.40890 


-0.29505 

0 

-0.50867 

0 

SJ, 

0 

-0.46861 

0 

0.27852 


0.059764 

0 

0.61154 

Q 

^6 

0 

0.098515 

0 

-0.61524 


0.001515 

0 

-0.18746 

0 

^8 

0 

-0.005411 

0 

0.25500 

"9 

0.000596 

0 

0.021697 

0 

^10 

0 

0.002858 

i 0 

i 

-0.028044 


The coefficients s^^q for the second and fourth modes were obtained by 

forcing the slope of Z^, as given by equation (BI5), to be zero at \|f = 

(that is, at x = 0). A similar table can be easily calculated for the 
pinned-edge plate . 
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TABLE I.- MODE-SHAFE FACTORS, EIGENVALUES, FREQUENCY 
RATIOS, AND STRUCTURAL INTEGRALS FOR 
FIRST FOUR NORMAL MODES 


(a) Plnned-edge plate®' 


Mode, 

n 

Ni 

N2 

Kn 

o^/o^ 

\ui 

Bnn 

1 

1 

— 

It 

1 

0.5 

0.57f2 

2 

1 

— 

2i( 

4 

.5 

2.0n2 

3 

1 

— 

3« 

9 

.5 

4.5n2 

4 

1 

— 

4n 

16 

.3 

8.0fl2 

1 


^A = B = 0 (m 4 n) 

TTin ■ TTin ' I / 


(b) Clamped-edge plate^^*^ 


Mo^, 

n 

% 

V , 

Kn 


V 

®nn 

1 

-0.629699 

-0.98250 

4.730 

1 

0.396 

4.88 

2 

-.66260 

-1.000778 

7.855204625 

2.7566 

.44o 

21.2 

3 

.7109645 

-.99997 

10.99560784 

5.4 o 4 

.506 

49.8 

4 

-.66120074 

-1.00000145 

14.13716549 

8.933 

.432 

76.4 


® (m ^ n) 

= B^jL = ^.362 
Bgj^ = B]^2 “ “S • 360 
Other B^ = 0 {m 4 n) 

c 

The significant figures shown were iised to avoid small -difference 
errors in mode shapes. 
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Stiffness parameter, cu,c/U 


(a) First -mode boiindary. 


Figure 1)-.- First-mode and second -mode stability bovindaries from a fouu*- 
mode analysis for clamped-edge panels for various values of structixral 

Po 

damping coefficient g. M = 1.3} f = — = 0. 




Mass ratio, \/ft 
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0 .1 .2 .3 .4 .5 .6 .7 

Stiffness parameter, tU|C/U 


(b) Second -mode boxindary. 
Figure 4.- Concluded. 




ratio, I //It 
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Stiffness parameter, o),c/U 


(a) First-mode and second-mode boundaries for g = 0. 



(b) First-mode boxindary for various values of g. 


Figure 6.- Stability boTxndaries from a two-mode analysis for clamped-edge 


panels for various values of g. M=1.5>f 


= Oj ^ = 0 


and 1.0. 




Mass ratio, 1/^ 
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Stiffness parameter, aijC/U 

Figure 7 .- Stability boimdaries from a two-mode analysis for clamped -edge 

Po 

panels for various values of f. M = 1.5; g = 0; — =1.0. 
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Stiffness parameter, o>,c/b 


(b) Pinned-edge panel. 


Figure 8.- Stability boundaries from a two-mode analysis for 
and pinned-edge panels for various values of g. M = |/2j 


clamped-edge 
















Mass ratio, l/u. 
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Figure 9*- 


Stability boundaries from two-mode and four-mode analyses for 

^ Po 

clamped-edge panels. M = 1.5o» g = f = -^ = 0* 
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Figure 11.- Structural damping coefficient g as a function of panel 
thickness ratio t/c from a two-mode analysis for clamped -edge 

Pn 

aluminiim panels in sea-level air. f = -^ = 0. (Refer to dashed 
hyperbolas in figs. 2, 8 (a), and 9 .) 
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NACA - Langley Field, Va. 


Figure 12.- Minimum panel thickness required to prevent flutter as 
function of ffech number for aluminum panels in sea-level air. 



